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ABSTRACT 

The shocks produced in the intergalactic medium during large-scale structure 
formation accelerate a population of highly relativistic electrons which emit synchrotron 
radiation due to intergalactic magnetic fields. In a previous paper (Loeb & Waxman 
2000) we have shown that these electrons cool primarily by inverse-Compton scattering 
of the microwave background photons and can thereby produce the observed intensity 
and spectrum of the diffuse 7-ray background. Here we calculate the intensity and 
angular fluctuations of the radio synchrotron background that results from the same 
high-energy electrons, as well as the expected angular fluctuations in the 7-ray 
background. On angular scales smaller than a degree, the predicted fluctuations in 
the microwave background temperature are of order 40;uK((^b/0.01)(i//10GHz)"^, 
where is the magnetic fraction of the post-shock energy density. This foreground 
might have already dominated the anisotropy signal detected in existing low-frequency 
CMB experiments, and can be identified with confidence through multi-frequency 
observations. Detection of the synchrotron fluctuations would allow to determine the 
strength of the intergalactic magnetic field. We predict a strong correlation between 
high-resolution maps taken at low-frequency radio waves and at high-energy 7-rays. 
Young X-ray clusters may also appear as radio or 7-ray clusters. The detailed study of 
this correlation will become easily accessible with the future launch of GLAST. 



Subject headings: Cosmology: diffuse radiation 



1. Introduction 

More than a third of all X-ray clusters with luminosities ^ 10^^ erg possess diffuse 
radio halos (Giovannini, Tordi, &: Feretti 1999). Based on energy arguments and circumstantial 
evidence, these radio halos are believed to be caused by synchrotron emission from shock- 
accelerated electrons during the merger events of their host clusters (Harris et al. 1980; Tribble 
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1993; Feretti 2000; see Liang et al. 2000 for references to alternative, less successful models) . These 
highly-relativistic electrons cool primarily through inverse-Compton scattering off the microwave 
background. Since their cooling time is much shorter than the dynamical time of their host cluster, 
the radio emission is expected to last as long as the shock persists and continues to accelerate fresh 
electrons to relativistic energies. Intergalactic shocks also occur along the filaments and sheets 
that channel mass into the clusters. These structures, also traced by the galaxy distribution (e.g. 
Doroshkevich et al. 1996), are induced by gravity and form due to converging large-scale flows in 
the intergalactic medium. 

In a previous paper (Loeb k, Waxman 2000), we have shown that most of the diffuse 7-ray 
background (Sreekumar et al. 1998) might have been generated by the shocks resulting from the 
formation of large-scale structure in the intergalactic medium. These shocks produce a population 
of highly-relativistic electrons with a maximum Lorentz factor ^ 10^ that scatter a small fraction 
of the microwave background photons in the present-day universe up to 7-ray energies, thereby 
providing the 7-ray background. The predicted flux agrees with the observed diffuse background 
over more than four decades in photon energy, provided that the fraction of shock energy carried 
by relativistic electrons is ~ 0.05 (a value consistent with that inferred for supernovae shocks). 
The brightest shocks arc caused by the accretion onto (or mergers of) X-ray clusters. The same 
electrons that emit 7-rays by inverse-Compton scattering of microwave background photons, also 
produce synchrotron radiation in the radio band due to intergalactic magnetic fields. The existence 
of magnetic fields with an amplitude ;^ 0.1/iG, is inferred in cluster halos (Kim et al. 1989; 
Fusco-Femiano ct al. 1999; Rephaeli et al. 1999; Kaastra et al. 1999), and is also required for the 
Fermi acceleration of these electrons. The appearance of radio halos around young X-ray clusters 
is therefore a natural consequence and an important test of our model for the extragalactic 7-ray 
background. The combination of radio and 7-ray data can be used to calibrate and determine 
the strength of the intergalactic magnetic field. 

The intergalactic synchrotron emission contaminates the microwave background anisotropics 
- relic from the epoch of recombination, and as such needs to be considered in the design 
and analysis of anisotropy experiments at low frequencies. Previous estimates of synchrotron 
contamination of the cosmic anisotropics focused on Galactic emission, which occurs primarily at 
low Galactic latitudes and large angular scales (Tegmark et al. 2000). In this Letter, we calculate 
the intergalactic synchrotron contribution to the fluctuations in the radio sky as a function of 
frequency and angular scale. We assume that most of the emission originates from the virialization 
shocks around X-ray clusters, and use the Press- Schechter mass function to describe the abundance 
of such clusters as a function of redshift. Our simplified analytic calculation should be regarded 
only as a first estimate of the synchrotron fluctuations and could be refined with future numerical 
simulations (Keshet et al. 2000). Throughout the paper we adopt the popular set of cosmological 
parameters (Ostriker &: Steinhardt 1995) with VLm = 0.3, = 0.65, Vtb = 0.05, h = 0.7, (Tg = 0.9, 
and n = 1. In §2 we describe the method of the calculation. In §3 we discuss the results and their 
main implications. Finally, we summarize our conclusions in §4. 
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2. Method of Calculation 

2.1. Emission from a Single Halo 

Let us first consider the luminosity produced by inverse-Compton and synchrotron emission 
of electrons accelerated to high energies in the accretion shock around a halo of mass M at a 
redshift z. Assuming that the shock is strong (i.e. the gas is not pre-heated), the population of 
shock-accelerated electrons acquires a number versus energy distribution with a power law index, 
•p = d\a.Nej d\\\.E = —2 ( |Bell 1978 , Blandford & Ostriker 1978). The inverse-Compton luminosity 
is then given by 

21n(7max)"m At"Zp 2 

where = + ^b, IJ-^p is the average particle mass in the gas including electrons (/i = 0.65), 
T is the post shock temperature, 7max ~ 10'' is the maximum Lorentz factor of the accelerated 
electrons (Loeb & Waxman 2000), M is the mass accretion rate across the virialization shock, 
and is the fraction of thermal shock energy carried by relativistic electrons. The synchrotron 
luminosity is then given by 

uLT = ^^^iMA.Lf, (2) 

^tCMB 

where ub and uqmb are the energy densities of the post shock magnetic-field and the cosmic 
microwave backround (CMB), respectively. 

The velocity dispersion, a, of a halo can be related to its mass through the relation 
where 

H = Hoa-^'^g{a) (4) 

is the Hubble parameter, with a = (1 + z)^^ and g{a) = + ^ko^ + (1 ~ — ^a)^]^^^- 
The numerical constant was chosen so that, for an isothermal sphere of velocity-dispersion 
cr, the quantity M is the mass enclosed within a sphere of average density p = 200/3c, where 
Pc = SH'^/SttG is the critical density at the corresponding redshift. On dimensional grounds, 
the ensemble- averaged mass accretion rate onto a halo of velocity dispersion a can be written as 
(White 1994) 

M(M,z) = /aec^^^, (5) 

while the shock temperature and radius are given by 

kT{M,z)= fTpmpa\M,z), (6) 



(7) 
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The factors /t, face and fsh are dimensionless numbers of order unity, which depend only weakly 
on the background cosmology for the relevant range of 0.3 ^ < 1, and should be calibrated 
based on numerical simulations (see, e.g. Miniati et al. 2000). We have chosen the normalization 
in equation so that for an isothermal sphere density distribution, the post shock density is 
Psh = 50/?c/3, similar to the result for the self-similar accretion solution in an Einstein de-Sitter 
{^Im = 1) cosmology (Bertschinger 1985). Note, that in this case the mass enclosed within r < r^h 
is 2M. 

We assume that the magnetic field energy density is a fixed fraction of the post-shock energy 
density, 

u,(M,.)=(,lkT(M.^)^^£l—. (8) 
With the above definitions, we find 



B{M,z) = 0.14 



f: 



n^J VO.Ol ' a'h%\-'[g{a)f^Mll'pG, (9) 

where Hq = 70/i7okm s"^ Mpc"\ M = lO^^MuM©. Thus, for values of Cb ~ 0.01 weh below 
unity, we obtain a magnetic field amplitude consistent with observations (Kim et al. 1989; 
Fusco-Femiano et al. 1999; Rephaeli et al. 1999; Kaastra et al. 1999) and with our model for the 
7-ray background (Loeb & Waxmann 2000). The implied synchrotron luminosity is 

uLT = 5.1 X lO^Vacc/^-^/Vsi;^ (^)' (^) (^) h^roa[9{a)? T^^ erg s^, (10) 

which can be related to the cluster mass through equation (^), 

Mi4 = 0.54/-'/'/ifo'a3/2[5(a)]-ir^/^, (11) 

where T^ev = (A;T/lkeV). While 7-ray observations can be used to calibrate ~ 0.05 (Loeb 
& Waxman 2000), radio synchrotron observations of clusters calibrate The synchrotron 



luminosity predicted by equation (10) is consistent with the values observed for the brightest radio 
halos of unrelaxed X-ray clusters for ~ 0.01 (see the upper envelope in Fig. 9 of Liang et al. 
2000, obtained for h = 0.5 and z ~ 0.3). Similarly to observations, our model predicts a radio 
luminosity which increases with cluster temperature and a nearly flat radio spectrum (see Fig. 6 
in Liang et al. 2000). 

The above expressions apply only to clusters which possess a strong shock with a Mach 
number T ^ 1. Non-accreting clusters or low-mass clusters which accrete gas that was already 
pre- heated close to their virial temperature (e.g., due to prior collapse of their environment or 
because of supernova energy injection), will not show significant non-thermal emission. For a shock 
front moving at a moderate Mach number, the accelerated electrons acquire a number-energy 
distribution with a power-law index p = {r + 2)/{r — 1), where r = 8T^/(6 -|- 2T^) is the shock 
compression ratio for a gas with an adiabatic index 7 = 5/3. This in turn yields steep non-thermal 
spectra with i/LJ^'^^^" oc z^^P/(^^~^)] ; implying that for weak shocks with T ^ 2, most of the 
non-thermal emission is at undetectable low frequencies. 
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2.2. Background Flux and Anisotropy 

The observed intensity at Earth of inverse-Compton and synchrotron emission is given by 
integrating (Q) and (^) over halo mass and cosmic time, t. Thus, 

where angular brackets denote sky average and dn/dM is the differential number density of halos 
per physical volume as a function of their mass. Using the Press-Schechter (1974) mass function, we 
find for the cosmological parameters we have chosen, = 1.5/acc/r(Ce/0-05)keV cm~^s~^sr~^. 

Thus, our model reproduces the 7-ray background intensity for a reasonable choice of parameters, 
namely /acc/T(Ce/0.05) 1 (Loeb & Waxman 2000). Given the calibration of the 7-ray background 
we predict a radio background with 

= 5.0 X lO-^Vacc erg cm-2 ,,-\ (13) 



/sh/ vo.osy vo.oi 

We find that 90% (50%) of the contribution to this background comes from massive clusters with 
M > 10^^ Mq (Af > 5 X lO^^M©), which are less sensitive to pre- heating of the accreting gas and 
for which the strong shock assumption is likely to apply. 

The intensity fluctuations can be calculated from 

52/(0) EE (4(O)/,(0))-(4)2 



d.^ /dM^ [^-^^(M,^)f 



dz J ^^'^ dM 167r3(l + z)^r^^{M, z) ^'^ 



edAiz) 



rsh{M, z) 



(14) 



Here, {Iy(fi)Iu{6)) is the sky average of the product of intensities along two lines of sight, 1 and 
2, separated by an angle 0, dA is the angular diameter distance, and Pi^2 is the probability that 
line-of-sight 2 passes through a halo of radius r^-^ given that line-of-sight 1 passes through the 
same halo, 

Pi,2{x) = - dcj) dyyxl' ^' ^ ' (15) 

vr Jo JO 1 0, otherwise . 



We have neglected in equation (14) halo-halo correlations, and assumed that the only contribution 
to the integral is from both lines-of-sights passing through the same halo. This assumption is 
justified since the optical depth (the average number of halos along a given line of sight) is much 
smaller than unity for the typical parameter values. In particular, we find that clusters with 
M > 5 X lO^'^M© (M > IO^^Mq) cover < 10% (< 2%) of the sky Note that the fractional 
fluctuations, ^^//(I)^ are independent of /t, /acc, Ce, and ^s, and their dependence on /sh is 

]/sh = l- 



simply given by f^S^ I {h^O) / {I)'' = [5^ I (9) / (I)^^ 
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3. Results 

The fluctuations in the radio and 7-ray background intensities, derived from equation (|l^, are 
shown in Figure 1. Although the synchrotron background amounts to only a small fraction of the 
CMB intensity, (/J")//^^^^ ^ g ^ io-6(/^//2j(^^/o.oi)(i./10GHz)-3 for /acc/rl^e/O.OS) = 1, 
its fluctuations could dominate over the primordial CMB fluctuations at low photon frequencies, 
^ 10 GHz. Our results imply that radio emission from cluster shocks contributes a fluctuation 
amplitude between ~ 30 and 50/iK xifr/ f^^){^B /O -01) {i^/WGRz)-^ to the CMB on angular 
scales between 1 and 0.1°, respectively. Interestingly, current anisotropy experiments are just 
sensitive to this level of fluctuation^. Existing detections by CAT (50 it 15/iK at 15 GHz on 
0.2-0.5°scales) and OVRO (56t|^/iK at 20 GHz on 0.1-0.6°scales), as well as 95% upper hmits 

40;uK on arcminute scales at 9-15 GHz by the ATCA, RYLE and VLA detectors) are consistent 
with our prediction^. 

For bright clusters, the flux increment due to synchrotron emission by the non-thermal 
electrons accelerated in the cluster accretion shock is larger at low frequencies than the Sunyaev- 
Zeldovich decrement due to the thermal intra-cluster electrons. For a gas density profile of an 
isothermal sphere, 

^ ^ -m..M^ {^) {^) (4) C'»-b,.r/»M- (^) 



lOGHz 

(16) 

In the X-ray regime, the central surface brightness due to thermal bremsstrahlung emission by the 
intra-cluster gas exceeds by orders of magnitude the inverse-Compton brightness; however, due to 
the rapid decline of the thermal brightness with projected radius outside the cluster core (oc 0~^), 
we find that the non-thermal brightness might dominate around the shock radius. 

The cumulative (all sky) number of radio and 7-ray halos as a function of flux or surface- 
brightness thresholds are shown at the top and bottom panels of Figure 2. We find that ~ 30 
7-ray halos have a flux exceeding the EGRET detection threshold of ~ 10~^^erg cm~^ s~^. Such 
halos may therefore constitute, as recently pointed out by Totani & Kitayama (2000), a significant 
fraction of the unidentified extra-Galactic EGRET sources (~ 60 sources over all sky; see Ozel &; 
Thompson 1996). Note, however, that since the angular extension of the brightest halos is large, 
a careful analysis is required to determine the detectability of such halos by the source search 
analysis which was applied to the EGRET data. 



'See table summary of current experiments at http://www.hep.upenn.edu/~max/index.html, or at http://cfa- 



www.harvard.edu/~mwhite/cmbexptlist.htm] 



* Note that although the cluster contribution declines rapidly at high frequencies, some of these experiments, like 
CAT or OVRO, measured fluctuations at only one frequency and could not reject a synchrotron contribution to the 
measured signal. 
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4. Conclusions 

The production of a fluctuating synchrotron background by strong intergalactic shocks is a 
natural consequence of structure formation in the Universe. The brightest emission originates from 
the virialization shocks on Mpc scales around newly formed, massive X-ray clusters. We estimate 
that the combined emission from all clusters produces fluctuations in the microwave background 
temperature of order 40//K(^5/0.01)(^e/0.05)(i//10GHz)~^ on sub-degree scales (Fig. 1), where 

~ 0.01 is the value required to explain the brightest radio halos of nearby clusters and ~ 0.05 
accounts for the 7-ray background. The foreground synchrotron fluctuations might be comparable 
to the anisotropy signals detected by existing low-frequency CMB experiments, and can be easily 
isolated through multi-frequency observations. Polarization anisotropy experiments could then 
constrain the coherence length of the intergalactic magnetic field. 

Our model predicts a fluctuation amplitude ^ 40% in the 7-ray background intensity on 
sub-degree scale (Fig. 1), and the existence of extended, ^ 1°, 7-ray halos, associated with 
newly formed massive clusters (Fig. 2). On scales larger than a degree the fluctuation amplitude 
declines and is well below the anisotropy limits from EGRET (see Fig. 5 in Sreekumar et al. 
1998), although somewhat higher than expected from an analogy to the X-ray background (Loeb 
Waxman 2000). Detection of the predicted signals will provide a calibration of the uncertain 
model parameter ^g. The high-energy maps required to detect the predicted anisotropy signal will 
be made between 20 MeV and 300 GeV by the GLAST instrument^ (planned for launch in 2005), 
which is expected to be more sensitive than EGRET by an order-of-magnitude ( [Bloom 1996[ ). 
The predicted 7-ray halos may constitute a significant fraction of the unidentified extra-Galactic 
EGRET sources (see also Totani & Kitayama 2000). However, since the angular extension of the 
brightest halos is large, a more careful analysis is required to assess the detectability of such halos 
by EGRET. 

A future, dedicated, all-sky anisotropy experiment, operating at several frequencies below 10 
GHz, would be able to map the fluctuations in the intergalactic synchrotron background. The 
resulting synchrotron map could then be cross-correlated with full-sky maps at hard X-ray or 
7-ray energies to confirm its cosmic origin. Identification of the synchrotron fluctuations together 
with their counterpart inverse-Compton emission of hard X-rays or 7-rays by the same population 
of shock-accelerated electrons, can be used to empirically determine the strength and spatial 
distribution of the intergalactic magnetic fields. Similarly, the correlation between radio and 7-ray 
halos may be detectable around individual X-ray clusters. Strong radio halos could be the best 
indicators for bright 7-ray clusters, which would provide the first obvious targets for GLAST. 

For the sake of simplicity, our model associated the intergalactic shocks with smooth spherical 
accretion of gas onto clusters, while in reality they result from asymmetric mergers as well as 
from converging flows in large scale sheets and filaments. The more realistic emission from these 



bttp:/ /glast. gsfc.nasa.gov/ 



for more details. 
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complex geometries can be best modeled through detailed hydrodynamic simulations (Miniati et 
al. 2000; Keshet et al. 2000). However, we note that mergers of comparable mass clusters would 
tend to produce only mild shocks due to the prior heating of the shocked gas, and hence result in 
negligible non-thermal emission due to the steep power-law slope of the accelerated electrons. 
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synchrotron (radio) background 




Fig. 1. — Fractional intensity fluctuations, ({Ii,{0)Ii^{9)) — in synchrotron (radio) 

and inverse-Compton (7-ray) background flux. The dimensionlcss coefficient, /sh, is of order 
unity [see definition in Eq. (7)]. The ratio of synchrotron to CMB intensity is = 
6 X 10~^(/r//s^h)(^B/0-01)(^/10Gmz)'~^, where fr is a dimensionless coefficient of order unity [see 
Eq. (6)], and the magnetic energy fraction is related to the magnetic field strength in Eq. (9). 
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Fig. 2. — The cumulative all-sky number of inverse-Compton (7-ray) and synchrotron (radio) 
halos with observed flux (top) or surface-brightness (bottom) exceeding a value of vfi, or vl,^, 
respectively. The normalization of the inverse-Compton halo curve is fixed by the requirement that 
the integrated emission would produce the observed 7-ray background. The dashed line represents 
the number of 7- ray halos with angular radius smaller than 1°. The synchrotron halo flux is 
proportional to and the curve shown in the flgure corresponds to = 0.01. The dimensionless 
coefficient, /sh, is of order unity [sec definition in Eq. (7)]; the index a = 2 for the inverse-Compton 
emission, and a = 4 for the synchrotron emission. 



